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ABSTRACT 
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initial  conditions,  may  be  determined  from  limited  measured  data.  The  use  of 
software  developed  to  remove  the  resulting  ship  motion  from  the  measured 
airwake  data  is  also  shown. 
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FFG-7  Ship  Motion  and  Airwake  Trial : 
Part  II  -  Removal  of  Ship  Motion  Effects  From 
Measured  Airwake  Data 


EXECUTIVE  SUMMARY 


A  helicopter/ship  computer  model  has  been  developed  to  simulate  the  complex  interactions 
in  the  dynamic  interface  between  ship  and  helicopter,  in  particular  between  the  ‘Adelaide’ 
class  FFG-7  frigate  and  Sikorsky  S-70B-2,  a  combination  used  by  the  Royal  Australian  Navy 
(RAN).  The  model  includes  helicopter  flight  dynamics  and  engine  dynamics,  undercarriage 
dynamics,  the  RAST  (Recovery  Assist,  Secure,  and  Traverse)  system,  ship  motion,  and  a 
representation  of  the  airwake.  The  original  ship  airwake  model  used  in  the  simulation  code 
has  been  studied  and  found  to  have  marked  differences  between  the  calculated  and  measured 
airwake  velocities,  when  applied  to  an  FFG-7  frigate. 


From  18  to  21  September  1989,  a  trial  was  undertaken  aboard  HMAS  Darwin,  an  ‘Adelaide’ 
class  frigate,  fitted  with  stabilisers  and  RAST  equipment.  The  objective  of  the  trial  was  to 
measure  the  ship  motion  and  airwake  over  the  flight  deck  for  a  variety  of  wind-over-deck 
velocities  for  use  as  a  data  base  in  the  helicopter/ship  simulation  code. 


Given  limited  measured  data  with  no  direct  measurement  of  ship  roll  and  pitch  attitude,  a 
procedure  has  been  developed  for  determining  ship  motion  and  removing  this  effect  from 
measured  airwake  data.  The  same  technique  for  determining  ship  motion  may  be  used  in 
future  trials,  with  a  notable  difference  being  that  if  attitude  is  measured  directly,  then  it  is  not 
necessary  to  use  derived  mean  values  as  an  approximation.  The  accuracy  of  the  calculated 
ship  velocity  may  be  improved  by  lengthening  the  data  recording  time  and  by  including  extra 
velocity  and  displacement  information  from  a  satellite  navigation  system. 


It  has  been  shown  that  the  mean  of  the  measured  airwake  is  approximately  the  same  as  the 
mean  of  the  airwake  when  corrected  for  motion  effects.  With  hindsight,  this  is  an  inevitable 
result,  as  the  mean  angular  rates  must  approach  zero  and  the  ship  is  assumed  to  have  mean 
linear  velocities  approaching  zero  in  the  motion  analysis.  For  the  applications  of  using  the 
measured  data  as  a  mean  flow  data  base  for  the  helicopter/ship  model,  or  for  comparison  with 
wind  tunnel  results,  the  procedure  for  removing  motion  effects  from  the  airwake  is  therefore 
not  warranted.  However,  it  should  be  noted  that  the  technique  for  determining  ship  motion 
developed  here  has  also  been  applied  to  determining  the  motion  of  a  ‘Perth’  class  DDG 
destroyer  and  for  calculating  the  position  of  the  flight  deck  of  HMAS  Jervis  Bay  during  a 
Black  Hawk  First  of  Class  Flight  Trial. 
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1.  INTRODUCTION 

In  1987,  the  Defence  Science  and  Technology  Organisation  (DSTO)  obtained  a 
helicopter/ship  computer  simulation  model  from  the  US  Naval  Air  Warfare  Center  (NAWC) 
Aircraft  Division  (previously  known  as  the  Naval  Air  Test  Center).  The  code  was  obtained 
under  the  auspices  of  The  Technical  Cooperation  Program  (TTCP).  The  model  can  be  used 
to  simulate  the  complex  interactions  in  the  dynamic  interface  between  ship  and  helicopter,  in 
particular  between  the  ‘Adelaide’  class  FFG-7  frigate  and  Sikorsky  S-70B-2,  a  combination 
used  by  the  Royal  Australian  Navy  (RAN).  The  model  includes  helicopter  flight  dynamics 
and  engine  dynamics,  undercarriage  dynamics,  the  RAST  (Recovery  Assist,  Secure,  and 
Traverse)  system,  ship  motion,  and  a  representation  of  the  airwake  (Ref.  1).  The  model  has 
been  substantially  modified  from  its  original  state  and  already  has  been  used  to  assess  a 
number  of  operational  problems  for  the  RAN  (Ref.  2).  Its  main  use  is  expected  to  be 
assisting  the  RAN  with  developing  safe  operating  envelopes  for  helicopters  operating  from 
ships  (Ref.  3). 

The  original  ship  airwake  model  used  in  the  simulation  code  has  been  studied  by  Erm 
(Ref.  4)  and  found  to  have  marked  differences  between  the  calculated  and  measured  airwake 
velocities  when  applied  to  an  FFG-7  frigate.  These  differences  are  attributed  to  the  origins 
and  subsequent  development  of  the  model,  in  which  measurements  based  on  different  class 
ships  were  used,  with  arbitrary  changes  then  made  to  improve  the  accuracy  of  a  real-time 
simulator. 

From  18  to  21  September  1989,  a  trial  was  undertaken  aboard  HMAS  Darwin,  an  ‘Adelaide’ 
class  frigate,  fitted  with  stabilisers  and  RAST  equipment.  The  objective  of  the  trial  was  to 
measure  the  ship  motion  and  airwake  over  the  flight  deck  for  a  variety  of  wind-over-deck 
velocities  (the  trials  instructions  are  detailed  in  Ref.  5)  for  use  as  a  data  base  in  the  acquired 
simulation  code.  The  instrumentation  package,  including  motion  sensors  and  anemometers, 
is  reported  in  Ref.  6  and  the  data  were  recorded  using  a  data  acquisition  system  (Ref.  7) 
developed  by  Air  Operations  Division  (AOD).  Prior  to  the  trial,  calibrations  were  obtained 
for  (a)  anemometers  in  the  AOD  low-speed  wind  tunnel,  (b)  angular  rates  on  a  rate  table,  and 
(c)  attitudes  on  a  tilt  table.  Just  before  sailing,  measurements  were  taken  while  the  ship  was 
alongside.  The  purpose  of  this  was  to  provide  procedural  practice  for  the  trials  team  and 
obtain  data  for  the  zero  ship  motion  case  as  a  baseline  when  investigating  effects  of  ship 
motion  for  the  ‘at  sea’  cases. 

The  results  of  the  trial  to  measure  ship  motion  and  airwake  are  documented  in  three  parts. 
Part  I  (Ref.  9)  gives  a  brief  description  of  the  scope  of  the  trial  and  details  of  the  data 
gathering  aboard  ship,  and  demonstrates  the  use  of  data  reduction  software  developed  to 
obtain  processed  data  in  imperial  units  from  the  raw  measured  data.  This  report  is  Part  II  and 
Part  III  (Ref.  10)  presents  the  airwake  mean  flow  and  turbulence  levels  over  the  flight  deck  of 
the  FFG-7. 

This  report  outlines  a  process  that  has  been  developed  whereby  the  ship  attitudes,  velocities, 
and  displacement,  including  initial  conditions,  may  be  determined  from  limited  data  through 
parameter  estimation.  The  use  of  software  developed  to  remove  the  calculated  ship  motion 
from  the  measured  airwake  data,  then  transform  it  from  ship  body  axes  to  ship-carried 
vertical  axes,  is  also  described.  Results  for  calculating  the  ship  motion  are  given  along  with  a 
comparison  of  airwake  data  in  ship  body  axes  with  airwake  data  in  ship-carried  vertical  axes 
that  have  the  motion-induced  component  removed. 

Imperial  units  are  adopted  in  this  report  because  (a)  they  are  used  exclusively  by  research 
workers  in  the  US  with  whom  AOD  is  collaborating  and  (b)  both  the  helicopter  and  ship 
referred  to  are  built  in  the  US  to  imperial  unit  specifications. 

2.  MEASUREMENTS  ACQUIRED  DURING  TRIAL 

The  airwake  over  the  flight  deck  was  measured  using  a  mobile  anemometer  mast  measuring 
over  30  feet  in  height,  that  was  placed  at  thirteen  positions  on  the  flight  deck  while  the  officer 
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of  the  watch  maintained  a  nominally  constant  relative  wind  over  deck  condition.  At  each  of 
the  thirteen  positions,  a  90  second  recording  run  was  made.  The  flight  deck  is  situated  at  the 
stem  of  the  FFG-7,  as  shown  in  Fig.  1,  and  the  mobile  airwake  mast  layout  and  positions  on 
the  flight  deck  are  shown  in  Fig.  2.  For  more  detailed  information  on  the  instmmentation  and 
methods  used  in  acquiring  the  data,  see  Ref.  9. 


Figure  1.  Flight  Deck  Location 
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Figure  2.  Mobile  Mast  Layout  and  Deck  Positions 


It  was  thought  that  the  height  of  the  mast  above  the  ship  centre  of  motion  was  such  that 
velocity  induced  by  the  ship  motion  would  be  significant.  This  induced  velocity  must  be 
removed  from  the  measured  data  and  the  resulting  airwake  velocity  data  transformed  from 
ship  body  axes  (as  measured)  to  ship-carried  vertical  axes  to  give  the  airwake  velocity  that 
would  be  seen  by  an  object  in  ship-carried  vertical  axes  over  the  flight  deck  (such  as  a 
helicopter).  Therefore  the  ship  linear  velocities,  angular  rates,  and  attitudes  must  be  known. 

The  quantities  measured  by  the  ship  motion  platform  included  the  three  acceleration 
components  (ax  ,  flJm,  azJ  and  roll,  pitch,  and  yaw  rates  (pm,  qm,  rm).  Two  additional 
quantities  recorded  from  the  ship  instrumentation  were  ship  heading  (yaw  attitude)  and  ship 
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speed.  It  should  be  noted  that  the  ship  speed  is  measured  using  a  device  which  would  tend  to 
average  the  actual  longitudinal  velocity  component.  As  reported  in  Ref.  9,  no  direct 
measurements  of  the  pitch  and  roll  attitudes  were  made,  although  some  indication  of  the 
average  attitudes  may  be  obtained  from  the  accelerometer  readings,  as  discussed  in  Section 
3.2. 

Thus  to  determine  the  ship  velocity  components  and  attitudes,  the  measurements  available 
are 

•  three  acceleration  components 

•  three  angular  rates 

•  average  longitudinal  velocity 

•  yaw  attitude 

Since  the  ship  was  not  at  a  known  trimmed  condition  when  recording  of  data  commenced, 
initial  conditions  are  unknown.  By  using  the  parameter  estimation  program  Compat 
(Section  3),  the  ship  attitudes,  velocities,  and  displacements,  including  the  initial  conditions, 
may  be  determined  after  first  making  a  few  reasonable  assumptions  as  discussed  in 
Section  3.2. 

The  parameter  estimation  program  Compat  may  also  be  used  to  check  measured  data  against 
expected  values  (as  determined  from  the  kinematic  equations  of  motion)  to  determine  bias 
(drift)  in  instrumentation. 

Data  reduction  of  the  measurements  obtained  from  the  trial  has  taken  place  in  two  distinct 
phases.  Phase  1  involved  deriving  the  ship  velocity  components  and  attitudes  from  the  ship 
motion  platform  instrumentation  using  program  Compat,  and  Phase  2  involved  correction  of 
the  anemometer  measurements  and  removal  of  the  ship  motion  components  from  the 
measured  airwake  velocities  using  program  AWSM. 

3.  CALCULATING  SHIP  MOTION 

To  determine  ship  motion  from  the  limited  data  measured  aboard  ship,  a  compatibility 
checking  technique  has  been  developed  using  the  program  Compat  (Refs  11  and  12). 
Because  of  the  intense  numerical  calculations  involved,  this  program  is  ran  on  a  RISC  6000 
machine.  The  technique  developed  here  has  also  been  applied  to  determining  the  motion  of  a 
‘Perth’  class  DDG  destroyer  (Ref.  13)  and  for  calculating  the  position  of  the  flight  deck  of 
the  training  ship  HMAS  Jervis  Bay  during  a  Black  Hawk  First  of  Class  Flight  Trial. 

The  axes  systems  used  in  the  motion  analysis  are  first  defined  in  Figure  3.  Note  that  at  time  t 
=  0,  the  origins  of  the  earth  axes  and  ship  body  axes  are  coincident,  but  the  initial  euler  roll 
and  pitch  angles  may  be  non-zero. 

Earth  Axes 

(origin  at  centre  of  ship  motion  at  time  t  =  0)  Ship  Body  Axes 


For  this  particular  trial,  program  Compat  is  first  used  to  determine  offset  errors  in  the 
instmmentation  then  calculate  motion  data,  namely  attitudes  and  velocities,  from  the  limited 
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measured  data.  This  section  describes  the  unique  subroutines  required  for  determining  ship 
motion,  and  the  procedure  to  be  followed  to  obtain  the  required  output. 

3.1  Program  Compat 

The  compatibility  checking  program  Compat  has  been  modified  from  its  original  form  to 
allow  it  to  run  on  personal  computers  (PCs)  with  limited  memory.  The  resulting  code  has 
been  rewritten  in  a  structured  way  so  that  the  code  is  much  easier  to  understand  and  modify, 
but  in  the  process  has  lost  the  ability  to  estimate  time-shifts  in  flight  data.1  However,  this 
version  of  the  code  is  more  easily  transported  to  different  machines. 

In  simple  terms,  Compat  takes  the  rigid  body  kinematic  equations  of  motion  in  differential 
form,  integrates  them,  and  then  compares  the  results  with  measured  data.  In  an  iterative 
process,  parameters  (that  in  this  case  are  instrumentation  offset  errors)  are  then  adjusted  to 
obtain  a  better  match.  If  the  measured  data  contained  no  offset  errors  then  there  would  be  no 
difference  between  the  measured  response  and  the  predicted  response  from  the  integrated 
equations  of  motion.  For  the  ship  motion  application  considered  here,  the  process  typically 
requires  25  iterations  until  convergence  occurs,  at  which  point  the  instrumentation  errors  and 
ship  motion  data  are  output.  The  optimisation  procedure  used  in  Compat  is  outlined  in 
Figure  4. 

For  the  ship  motion  application,  the  control  inputs  are  taken  to  be  the  measured  accelerations, 
angular  rates,  and  angular  accelerations.  The  experiment  consists  of  the  assumed  velocities 
and  attitudes  resulting  from  these  inputs  (see  Section  3.2),  and  the  mathematical  model  is  the 
kinematic  equations  of  motion  with  allowances  incorporated  for  instrumentation  offset  errors. 


Control  Measured 


Updated  Parameters 


Figure  4.  Parameter  Estimation  Procedure  Used  in  Program  Compat 

Compat  requires  two  data  files  as  input.  The  COMDAT .  xxx  file  controls  operation  of  the 
program  by  defining  which  parameters  are  to  be  determined,  how  many  controls,  states,  and 
measurements  are  involved,  and  a  number  of  other  parameters,  defined  in  Appendix  A.  For 
the  ship  motion  analysis,  a  two-stage  method  is  required  using  two  COMDAT .  xxx  files  in  two 
separate  runs  (see  Section  3.3),  which  are  also  given  in  Appendix  A.  The  second  data  file 
required  by  Compat  (xxxxxxxx_Out)  contains  a  time  history  of  the  control  inputs  and 
measured  response  data  (see  Section  3.2). 

Three  output  files  are  created  by  Compat.  They  are  COMOUT .  FFG,  THIST .  DAT  and 
var .  dat.  File  COMOUT .  FFG  contains  run-time  information,  most  details  of  which  will  not 
be  discussed  here.  File  THIST .  DAT  contains  the  measured  and  observed  responses 
(calculated  in  Resp .  f ,  see  below)  in  the  form  of  time  histories.  File  var .  dat  contains 
the  measured  ship  longitudinal  velocity,  the  calculated  states,  and  the  angular  rates  corrected 
for  instrument  bias  errors.  Figure  5  shows  the  program  modules  and  input  and  output  files 


1  The  time-shift  capability  was  removed  because  the  way  it  was  originally  implemented  was  too  memory 
intensive  for  a  PC.  It  is  intended  to  reintroduce  the  time-shift  capability  at  a  later  date. 
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and  Table  1  summarises  the  contents  of  each  data  file  and  module  that  is  discussed  further 
below. 

The  program  itself  has  three  problem-specific  subroutines.  These  subroutines  are  kept  as 
separate  modules  (called  Deriv .  f ,  Init .  f ,  and  Resp .  f )  and  are  given  in  Appendix 
B,  along  with  the  executable  file  link,  which  compiles  and  links  the  required  modules  to 
create  Compat.  Subroutine  Deriv .  f  contains  the  derivative  equations  of  motion  and  the 
instrumentation  error  models,  Init .  f  contains  the  initial  conditions,  and  Resp .  f  contains 
the  calculated  response.  The  contents  of  these  subroutines  are  discussed  further  below. 

For  the  ship  motion,  the  control  inputs  (Fig.  4)  used  are  the  angular  rates  and  accelerations 
and  the  linear  accelerations.  An  error  model  for  the  rates  and  linear  accelerations  is  also 
included,  which  allows  an  offset  in  the  measured  data.  Equations  below  are  in  file  Deriv .  f 
(Appendix  B). 

Input  Files  Output  Files 


Problem  Specific 


Figure  5.  Program  Modules  and  Input  and  Output  Files 


Table  1.  Summary  of  Compat  Modules  and  Data  Files 


Module/File 

Description 

COMDAT . xxx 

Program  operation  inputs 

xxxxxxxx_Ou  t 

Time  history  of  control  inputs  and  measured  response  data 

COMOUT . FFG 

Run-time  information,  including  estimates  of  parameters 

THIST.DAT 

Time  history  of  measured  data  and  calculated  responses 

var . dat 

Time  history  of  calculated  (true)  ship  motion 

Compat  (Executive) 

Executive  routine  controlling  execution  of  program 

Comsub 

Matrix  manipulation  routines 

Rk4new 

Integration  routine 

Init 

Determines  initial  conditions 

Deriv 

Derivative  equations  and  error  models  (math  model  in  Fig.  4) 

Resp 

Predicted  responses  (calculated  ship  motion) 

The  calculated  roll,  pitch,  and  yaw  rates,  p,  q,  and  r,  are  given  by 

P  =  Pm+Pb 
*7  —  Qm  "b  Qb 
r  =  rm  +  rb 
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where  control  inputs  pm,qm,rm are  the  measured  rates  and  parameters  pb,qb,rbaxo  the 
instrument  bias  (offset)  errors  in  measured  rates. 

The  calculated  linear  accelerations  along  the  x,  y,  and  z  axes,  ax,ay,az  ,  are  given  by 

ax  =  ax  +  axb 
ay=ay+  dyb 

az~az  +  azb 

where  parameters  aXb,ayb,aZb  are  the  offset  errors  in  measured  linear  accelerations  and 
a'x,ay,a'z  are  the  measured  accelerations  corrected  for  accelerometer  offset  from  the  centre  of 
ship  motion  and  gravitational  effects  (Ref.  14),  given  by 

a'x  =  axm  +  (<?2  +  r2)x  -  ( pq  -  r)y  -  ( pr  +  q)z~g  sin  9 

a'y  =  aym  -(pq  +  r)x  +  (p2 +  r2)y-(qr- p)z  +  g cos  9 sirup 
a'=a7  -  ( pr  -  q)x  -  ( qr  +  p)y  +  (p2  +  q2)z  +  geos  9c os  <p 
where  control  inputs  aY  ,av  ,a7  are  the  measured  linear  accelerations,  p,q,r  are  the 

r  Am  Jm  4m 

calculated  roll,  pitch,  and  yaw  accelerations,1  and  parameters  x,y,z  are  the  linear 
displacements  of  the  accelerometers  from  the  centre  of  ship  motion. 

The  state  variables  (referred  to  as  states)  defining  the  system  are  the  Euler  angles  ( <p,9,yr), 
linear  velocities  in  ship  body  axes  ( U,V,W ),  and  displacements  in  earth  axes  (X,Y,Z), 
which  are  found  by  integrating  the  following  standard  rigid  body  kinematic  equations  of 
motion  (Ref.  14): 

0  =  p  +  qsiruptan9  +  rcos(ptan9 

9  =  qcos<p-rsirup 

yr  =  (q  sin  (p  +  r  cos  (p)  sec  9 

U  =  ax-  qW  +  rV 

V  =  ay-rU  +  pW 
W  -  az-  pV  +  qU 

X  =  Ucos9cosy/+  V(sin  (p  sin  9  cos  if/  -  cos  (p  sin  yr)  +  W(cos  (p  sin  9  cos  yf  +  sin  (p  sin  y/) 

Y  =  Ucos9siny/  +  E(sin  (p  sin  9  sin  y/ + cos  (p  cos  y/)  +  W(cos<p  sin  9  sin  y/  -  sin  (p  cosy/) 
Z--Usir\9  +  V  cos  9  sin  (p  +  W  cos  9  cos  (p 

When  the  data  runs  were  recorded,  the  initial  conditions  for  many  of  the  states  were 
unknown.  File  Init .  f  (Appendix  B)  gives  the  user  several  options  for  determining  initial 
conditions.  These  include  a  least  squares  fit  to  the  first  few  data  points,  updating  the 
measured  value  using  error  models  defined  in  Resp .  f ,  or  estimating  the  initial  condition  as 
separate  parameters.  Since  U,  V,  W,  Z,  <p,  and  9  were  not  measured,  the  only  option  is  to 
estimate  initial  conditions.  By  definition,  the  initial  conditions  for  A,  Y,  and  y/are  zero. 

The  observed  responses  which  are  compared  with  measured  data  are  defined  in  subroutine 
Resp  .  f  (Appendix  B)  and  are  generally  the  integrated  values  of  the  state  derivatives,  with 
the  exception  of  X  and  Y.  Initially,  the  measured  data  for  X  were  assumed  to  be  Ut ,  and  for  Y 
were  assumed  to  have  a  mean  value  of  zero  since  no  measured  data  were  available  for 
comparison.  However,  this  resulted  in  parameters  which  were  obviously  incorrect.  A 
mechanism  was  required  whereby  the  derivative  equations  for  X  and  Y  could  still  be  used  to 


1  The  true  angular  accelerations,  p,q, and  r ,  are  input  to  the  program  as  control  inputs  UU(7),  UU(8),  and 
UU(9)  (Appendix  B).  They  are  derived  from  the  filtered  measured  rates  (Ref.  9),  which  are  assumed  to 
have  only  an  offset  error,  and  thus  are  assumed  to  be  true  accelerations. 
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calculate  X  and  Y  while  avoiding  optimising  the  parameters  to  match  measured  and  observed 
responses  for  X  and  Y.  A  solution  was  to  set  the  observed  responses  for  X  and  Y  equal  to  the 
measured  data,  which  were  arbitrarily  set  to  zero. 

3.2  Preparing  Time  History  Input  Data  for  Program  Compat 
The  time  history  data  file  (xxxxxxxx_Out),  which  is  used  as  input  to  Compat,  must  be  in 
the  format  defined  in  Table  2. 

The  parameters  measured  by  the  ship  motion  platform  (Ref.  9)  included  the  three  acceleration 
components  (ax  ,ay  ,az  )  and  roll,  pitch,  and  yaw  rates  ( pm,qm,rm ).  Two  additional 

parameters  recorded  from  the  ship  instrumentation  were  ship  heading  (y/m)  and  ship  speed 
( Um ).  It  should  be  noted  that  the  ship  speed  is  measured  using  a  device  which  would  tend  to 
average  the  actual  longitudinal  velocity  component. 


Table  2.  Input  Format  Required  by  Compat 


Column 

Variable 

Comments 

1 

2-4 

Time 

ax  ,  ,  ciy 

jm  -wn 

5-7 

Pm '  Qm  >  rm 

•  Control  Inputs 

8-10 

p,q,r 

1 

11-13 

U  V  w 

v rrv> rvm 

I 

14-16 

>  Measured  Responses 

17-19 

X  Y  7 

1  m^m 

As  reported  in  Ref.  1,  no  direct  measurements  of  the  pitch  and  roll  attitudes  were  made, 
although  some  indication  of  the  mean  attitudes  may  be  obtained  from  the  accelerometer 
readings.  For  example,  in  a  static  condition,  if  the  ship  rolled  90°,  the  lateral  accelerometer 
would  read  lg  (close  to  32.2  ft/s2).  Thus  attitude  derived  from  accelerometer  measurements 
is  given  by 


( pm  —  -  arcsin 


ym 


0m  -  arcsin 


V  s  ) 


However,  for  the  dynamic  case  of  ship  motion,  the  accelerometers  also  sense  linear 
accelerations  of  the  ship  due  to  wave  movement.  Since  the  ship  is  maintaining  a  steady 
forward  velocity,  it  may  be  assumed  that  the  mean  acceleration  components  due  to  wave 
movement  are  zero.1  Thus  the  mean  values  derived  from  the  accelerometers  correspond  to 
the  mean  attitudes,  as  determined  by  Compat,  and  so  provide  additional  information 
necessary  for  the  parameter  estimation  process  to  succeed. 

Lateral  and  vertical  velocities  in  body  axes  (Vm  and  Wm  respectively),  and  displacement  in 
earth  axes  ( Zm  ),  were  not  measured,  but  the  data  file  for  Compat  must  contain  values  to 
compare  with.  If  the  ship  is  sailing  at  a  constant  speed  and  heading,  then  it  is  reasonable  to 
assume  that  the  mean  velocities  and  vertical  displacement  over  a  sufficiently  long  time  period 
(taken  as  over  60  s,  see  Section  5.1)  will  be  effectively  zero.  Thus  by  using  constant  values 
of  zero  for  Vm  ,  Wm  ,  and  Zm  as  input,  Compat  will  calculate  F,  W,  and  Z  such  that  the  mean 
will  be  zero.  Constant  values  for  Xm  and  Ym  are  also  used  as  described  in  Section  3.1. 


1  During  this  trial  the  sea  conditions  were  relatively  benign.  If  the  sea  was  such  that  the  ship  forward  velocity 
was  not  reasonably  steady  during  the  data  recording  period,  then  either  the  recording  period  would  need  to 
be  increased  or  another  technique  would  need  to  be  developed. 
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In  order  to  create  the  xxxxxxxx_Out  input  file,  which  contains  both  measured  data  and 
data  derived  using  the  above  assumptions,  several  utility  programs  (written  in  Fortran)  are 
available  for  use  on  Macintosh  computers. 

Program  MacTrans  (Ref.  8)  is  first  used  to  obtain  data  in  imperial  units  from  the  data  file 
generated  using  program  MacShipRefine  (Ref.  9).  Appendix  C  includes  a  copy  of  file 
MOTION .  BLK,  which  contains  a  list  of  the  channels  required  in  the  correct  order  (LI ,  L2)  for 
use  with  other  utility  programs  outlined  below.  Before  running  MacTrans,  file 
MOTION .  BLK  should  be  renamed  to  TRANS .  BLK.  An  example  demonstrating  the  use  of 
MacTrans  is  shown  below  for  the  data  file  2  0  091158.DAT  (bold  italic  type  being  user 
inputs  and  <CR>  denoting  a  carriage-return). 

[TRANS  version  date  12-FEB-91] 

GO  STRAIGHT  TO  EVENT  LOOP? 

I/P  FILENAME  =  20091158 

Wind  60  deg,  20  kn  -  Position  10 

File  20091158  (filtered) 

I/P  FILE  RECORDED  ON  17-FEB-93  AT  10:11:28 

INTEGN  INT  =  0.0000E+00;  RUN  CPU  TIME  =  3  MIN.  24.47  SEC. 

TIME  FROM  0 . 0000E+00  TO  8.9500E+01  IN  STEPS  OF  5.0000E-02 

*  PRC 

PRINTING  IN  COLUMNS  : 

BLKS 

L1,L2 

<CR> 

IS  O/P  TO  TTY  REQRD  :  N 

*GOE 

**  RUNNING  ** 

*EXI 

The  resulting  data  file,  20091 158.COL,  is  then  used  as  input  to  program 
FFGConvert  (Attitudes).  This  utility  program  converts  the  MacTrans  column  format  into  the 
required  multi-column  form  (Table  2),  which  is  also  more  easily  read  in  by  other  programs. 
When  running  program  FFGConvert  (Attitudes),  a  window  appears  prompting  for  the  name 
of  the  MacTrans  output  file  (Fig.  6).  After  selecting  the  appropriate  file,  which  for  this 
example  is  20091 158.COL,  the  output  file  that  is  automatically  named  after  the  input  file, 
20091158  (Angles) ,  is  produced. 

Program  KaleidaGraph1  is  used  to  read  in  data  file  20091158  (Angles)  and  then  to  plot 
columns  14  and  15  ( <pm  and  0m).  The  mean  attitudes  may  then  be  determined  using  a  linear 
curve-fit.  The  equation  for  a  linear  curve-fit  is  used,  rather  than  a  constant  mean  value,  to 
allow  for  a  drift  in  the  mean  over  time.  The  results  for  roll  and  pitch  ‘attitudes’  are  shown  in 
Figure  7. 


1  ©  Abelbeck  Software,  distributed  by  Synergy  Software  (PCS  Inc.),  2457  Perkiomen  Avenue,  Reading,  PA 

19606,  USA. 
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Program  FFGConvert(heave)  generates  the  data  file  20091158  Out,  which  contains  zero 
values  for  the  lateral  and  vertical  velocity  components  and  vertical  displacement,  and  also 
calculates  values  of  0  and  6  according  to  the  linear  curve-fits  given  by  KaleidaGraph.  When 
running  FFGConvert(heave),  a  window  appears  (Figure  8)  which  prompts  for  the  offsets  and 
gradients  as  determined  above. 


=n -  •  === — output  from  FFGConuert(heauu)  =========== 

=P1=j 

Input  file  <MacTrans  ".COL") 

Input  Attitude  Equation  Coefficients  <deg> 

Roll  Offset:  -.67813 

Roll  Gradient:  -.0026244  Y 

Pitch  Offset:  -.33838  1 

Pitch  Gradient:  .00012553 

STOP 

Figure  8.  Window  For  Running  FFGConvert(heave) 


3.3  Running  Compat 

Having  created  the  Compat  input  file  2  0  091158  Out  on  a  Macintosh  computer  (Section 
3.2),  the  file  is  transferred  to  a  Rise  6000  machine  using  ftp  in  the  text  mode  via  program 
Versaterm  Pro  (©  Abelbeck  Software).  On  the  Unix  based  machines,  spaces  in  filenames  are 
converted  to  underscores;  thus  file  20091158  Out  becomes20091158_0ut. 

It  was  determined  at  an  early  stage  that  Compat  would  not  always  converge  if  an  attempt  was 
made  to  optimise  simultaneously  all  the  parameters  required  for  all  the  states.  Because  the 
equations  for  the  Euler  angles  are  decoupled  from  the  other  equations,  the  parameters 
associated  with  the  angles  are  first  determined  (offset  errors  in  rates  and  initial  conditions  of 
attitudes).  The  parameters  are  then  set  constant  for  a  second  run  of  Compat  when  the 
parameters  associated  with  the  other  state  variables  are  found.  The  parameters  are  referred  to 
by  the  designated  numbers  defined  in  the  input  files  given  in  Appendix  A.  The  computer 
variable  names  are  often  slightly  different  from  the  nomenclature  defined  in  the  equations 
(Section  3.1).  Table  3  defines  the  relevant  parameters,  which  are  discussed  further  below. 

COMDAT  .0  02  (Appendix  A)  is  used  as  input  for  the  first  run,  where  parameters  4,  5,  6 
( Pb’%’rb  in  rad/s)  and  28  and  29  ( $0  and  Oo  in  rad)  are  estimated.  Parameters  31,  32,  33 
(the  position  x,y,z  in  ft  of  the  accelerometers  in  ship  body  axes)  are  set  to  constant  values. 
A  mechanism  to  ensure  that  the  motion  parameters  are  not  optimised  using  the  equations  for 
X,  Y,  and  Z  (subroutine  Derivs,  Appendix  B)  was  to  use  weighting  parameters  42,  43,  and 
44  which  are  set  to  zero.  The  state  weighting  matrix  (Dl)  is  also  set  with  large  numbers  for 
the  Euler  angles  and  small  ones  for  the  other  states.  This  has  the  effect  of  conditioning  the 
parameters  so  that  they  will  be  weighted  towards  optimisation  by  only  using  the  angular  state 
equations.  An  example  of  running  Compat  is  shown  below. 


<blackjack  /arneya/compat>  Compat 
COMDAT. xxx  File  Number? 

2 

COMDAT .002  opened . 

Input  File? 

20 09115 8_Out 

20091158_0ut  opened. 

*****  nno  >  Data  file  size  -  Data  array  reduced  from  2500  to  1790 
points 

ITERATION  NUMBER  0 

ITERATION  NUMBER  1 
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ITERATION  NUMBER  2 

ITERATION  NUMBER  3 

ITERATION  NUMBER  23 

ITERATION  NUMBER  24 

STOP 


STOP 

<blackjack  /arneya/compat> 

Table  3.  Equation  Nomenclature  and  Equivalent  Computer  Symbols 


Parameter 

Number 

Equation 

Nomenclature 

Computer 

Symbol 

Comments 

l 

axb 

BAX 

2 

a>'b 

BAY 

►  Accelerometer  offsets 

3 

% 

BAZ 

4 

Pb 

BP 

1 

5 

% 

BQ 

•Angular  rate  offsets 

6 

rb 

BR 

- 

25 

U0 

UO 

26 

Vo 

VO 

•  Initial  velocities 

27 

Wo 

WO 

- 

28 

$o 

PHIO 

1 

■  Initial  attitudes 

29 

0o 

THO 

J 

31 

X 

BXCG 

32 

y 

BYCG 

•  Accelerometer  position  in  ship  body  axes 

33 

z 

BZCG 

■ 

38 

Zo 

Z0 

Initial  vertical  displacement 

42 

WX 

43 

- 

WY 

•  Weighting  parameters 

44 

- 

WZ 

- 

Note  that  the  number  of  data  lines  read  in  is  set  to  2500  (2500*0.05  =  125  s  data  run)  in 
COMDAT  .002,  but  Compat  will  automatically  adjust  for  data  files  smaller  than  this.  File 
COMOUT .  FFG  contains  run-time  information,  most  details  of  which  will  not  be  discussed 


11 


here.  The  part  of  COMOUT .  FFG  of  particular  interest  at  this  stage  is  towards  the  end  of  the 
file  where  the  estimated  parameters  are  given  along  with  the  Cramer-Rao  (C-R)  bound.  The 
relevant  part  of  the  file  is  duplicated  below. 


PARAMETER  STATUS 

A. P.  WEIGHT 

ANSWER 

C-R  BOUND 

4 

1 

.000000  99999.000 

-.020238 

.000012 

5 

1 

.000000  99999.000 

.000979 

. 000007 

6 

1 

.000000  99999.000 

. 000044 

. 000001 

28 

1 

.000000  99999.000 

-.006754 

.000638 

29 

1 

.000000  99999.000 

-.008442 

.000358 

As  a  general  rule,  if  the  C-R  bound  is  less  than  one  tenth  of  the  estimated  parameter 
(ANSWER),  then  a  high  level  of  confidence  can  be  attached  to  the  estimate.  As  can  be  seen 
above,  a  very  high  degree  of  confidence  is  associated  with  the  estimates  of  the  angular  rate 
offsets,  pb,  qb ,  and  rb,  (4,  5,  and  6),  although  the  confidence  level  for  the  initial  attitudes,  <f>0 
and  d0,  is  lower.  This  is  because  the  attitude  initial  conditions  are  obtained  from  mean  data 
rather  than  measured  data.  Figure  9  shows  a  comparison  of  the  measured  rates  with  the 
calculated  rates  determined  by  Compat,  and  Figure  10  shows  a  comparison  of  the  measured 
attitudes  (derived  from  accelerometers  for  roll  and  pitch,  compass  for  yaw)  with  the 
calculated  attitudes. 

COMDAT  .004  (Appendix  A)  is  used  for  the  second  run  where  parameters  1,  2,  and  3  (offsets 
in  accelerometer  measurements  in  ft/s^),  25,  26,  and  27  (initial  U,  V,  and  W  in  ft/s),  and  38 
(initial  Z  in  ft)  are  estimated.  Parameters  31,  32,  and  33  remain  constant  (as  for 
COMDAT .  002),  but  parameters  42,  43,  and  44  are  set  to  unity  so  that  the  derivative  equations 
for  X,  Y,  and  Z  are  used  (subroutine  DERIVS,  Appendix  B).  Parameters  4,  5,  6,  28,  and  29 
are  set  constant  to  the  values  determined  in  the  previous  run  using  COMDAT .  002  as  input. 
The  state  weighting  matrix  (Dl)  should  also  be  set  with  small  numbers  for  the  Euler  angles 
and  large  numbers  for  the  other  states.  An  example  of  running  Compat  is  shown  below. 


<blackjack  /arneya/compat>  Compat 
COMDAT. xxx  File  Number? 

4 

COMDAT .004  opened . 

Input  File? 

20091158_0ut 

20091158_0ut  opened. 

*****  nno  >  Data  file  size  -  Data  array  reduced  from  2500  to  1790 


points 

ITERATION  NUMBER  0 

ITERATION  NUMBER  1 

ITERATION  NUMBER  23 

ITERATION  NUMBER  24 

STOP 


<blackjack  /arneya/compat> 
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Again,  examining  file  COMOUT .  FFG  gives  the  parameters  and  the  Cramer-Rao  bounds. 


PARAMETER  STATUS 

A.P.  WEIGHT 

ANSWER 

C-R  BOUND 

1 

1 

.109218  99999.000 

-.055403 

.000386 

2 

1 

1.247285  99999.000 

.777414 

.000677 

3 

1 

-.120199  99999.000 

-.316788 

.000107 

25 

1 

7.829200  99999.000 

19.176580 

.020535 

26 

1 

.000000  99999.000 

.626917 

.035006 

27 

1 

.000000  99999.000 

-1.252886 

. 004925 

38 

1 

-2.317998  99999.000 

1.839305 

.096327 

It  can  be  seen  from  the  Cramer-Rao  bounds  that  all  the  estimated  parameters  have  a  high 
confidence  level.  The  final  product  (i.e.  the  ship  motion)  is  given  in  file  var .  dat ,  which  is 
now  transferred  back  to  a  Macintosh  using  ftp  in  the  text  mode  via  program  Versaterm  Pro, 
so  that  the  motion  effects  may  be  removed  from  the  measured  data. 

It  should  be  noted  at  this  stage  that  when  running  Compat  for  a  number  of  different  data  files 
(xxxxxxxx_Out  files),  any  COMDAT  .002  file  set  up  for  ship  motion  may  be  used,  but  the 
COMDAT  .004  file  must  be  edited  with  the  specific  parameters  4,  5,  6,  28,  and  29  found  in  a 
previous  run  using  COMDAT  .002. 
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Rate  Pitch  Rate  Roll  Rate 


Yaw  Attitude  Pitch  Attitude  Roll  Attitude 
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Figure  10.  Measured  Attitudes  Compared  With  Calculated  Attitudes 
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4.  REMOVING  SHIP  MOTION  FROM  AIRWAKE  DATA 

Ship  motion  effects  are  removed  from  the  measured  anemometer  data  using  a  program  on  the 
Macintosh  (written  in  Fortran)  called  AWSM  (Airwake  Without  Ship  Motion).  This  section 
describes  how  AWSM  removes  these  effects  and  gives  examples  of  running  AWSM  and 
associated  utility  programs. 

4.1  Program  A  WSM 

AWSM  reads  in  calculated  ship  motion  (from  Compat,  see  Section  3)  and  measured  airwake 
data  (from  MacShipRefine,  see  Ref.  9),  removes  the  motion  effects,  then  outputs  the 
corrected  airwake  data  to  file  xxxxxxxx  Output.  The  position  of  each  of  the 
anemometers  relative  to  the  centre  of  ship  motion  is  first  determined,  depending  on  the 
position  of  the  mobile  anemometer  mast  (Fig.  2).  Table  4  gives  the  displacement  of  the  base 
of  the  mobile  mast  for  each  of  the  mast  positions,  as  well  as  the  position  of  the  cup  and 
aerovane  anemometers. 


Table  4.  Mobile  Mast  Base,  Cup,  and  Aerovane  Anemometer  Positions 


Anemometer 

xbase  W 

Position  1 

-118.7 

-17.1 

-15.3 

Position  2 

-118.7 

-10.5 

-15.3 

Position  3 

-118.7 

0.0 

-15.3 

Position  4 

-118.7 

10.5 

-15.3 

Position  5 

-118.7 

17.1 

-15.3 

Position  6 

-142.5 

-17.1 

-15.8 

Position  7 

-142.5 

-10.5 

-15.8 

Position  8 

-142.5 

0.0 

-15.8 

Position  9 

-142.5 

10.5 

-15.8 

Position  10 

-142.5 

17.1 

-15.8 

Position  1 1 

-163.5 

-10.5 

-16.2 

Position  12 

-163.5 

0.0 

-16.2 

Position  13 

-163.5 

10.5 

-16.2 

Aerovane 

-181.9 

2.1 

-29.93 

Stbd  Cup 

-181.9 

21.7 

-21.7 

Port  Cup 

-181.9 

-21.7 

-21.7 

The  anemometers  on  the  mast  are  offset  from  the  base  along  the  X  and  Z  axes  as  shown  in 
Table  5.  The  position  of  each  of  the  anemometers  in  the  array  relative  to  the  ship  centre  of 
motion  is  calculated  as  follows: 


Table  5.  Offsets  in  Anemometers  Relative  to  Mast  Base 


Top 

Mid 

Low 

Hat  («) 

0.0 

1.48 

1.48 

xvert  (ft) 

0.0 

1.48 

1.48 

^Zortg(ft) 

-31.5 

-21.0 

-10.5 

Zta(ft) 

-31.83 

-21.33 

-10.83 

For  each  anemometer  array  (top,  mid,  and  low), 

^-Hat  ~  xbase  xlat 
A* vert  ~  xbase  xvert 
^Hlong  ~  ybase 
A yvert  ~  ybase 
A? long  ~  Zbase  Hong 
A^ lat  ~  Zbase  Hat 

The  corrected  anemometer  velocities  in  ship  body  axes  for  each  array  are  given  by 

uc  —  um  ~  ^Uship  +  Q^H-long  ~  r^yiong 

Vc  =  vm  -  A Vghip  +  rAxiat  -  pAzlat 
=  wm  ~  A Wship  +  pAyvert  -  qAxvert 

where 

A  U ship  —  U—  Umean 
Ay,.  =V-V 

^ v  ship  v  Y  mean 

A W  u-  =W-W 

Li,,ship  ''  rr  mean 

Note  that  the  incremental  ship  velocities  (fluctuations  about  the  mean  velocity),  rather  than 
the  actual  ship  velocities,  are  subtracted  from  the  measured  airwake  velocities.  This  is  to 
allow  for  the  case  where  all  the  measured  airwake  velocity  is  due  to  the  ship  velocity.  For 
example,  in  still  air  with  an  actual  ship  velocity  of  10  ft/s,  the  measured  airwake  will  be 
10  ft/s.  Thus,  if  the  actual  ship  velocity  was  subtracted,  the  corrected  airwake  velocity  would 
be  zero  which  is  clearly  incorrect.  It  should  also  be  noted  that  one  of  the  key  assumptions  in 
determining  the  ship  motion  using  Compat  is  that  the  mean  lateral  and  vertical  ship  velocities 
in  body  axes  are  zero. 

AWSM  also  transforms  the  velocities  from  ship  body  axes  to  ship-carried  vertical  axes,  which 
is  the  desired  format  for  comparison  purposes  with  wind  tunnel  measurements  and  for  use  as 
an  airwake  data  base  in  the  computer  simulation.  The  ship-carried  vertical  axes  origin  is 
coincident  with  the  ship  body  axes  (Fig.  1 1),  but  the  Z  axis  is  aligned  vertically  downwards, 
i.e.  along  the  local  g  vector. 


Figure  11.  Ship-carried  Vertical  Axes  (X’Y’Z’,  Z’  Parallel  to  Ze) 
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The  matrix  used  to  transform  from  ship  body  axes  to  ship-carried  vertical  axes  is  shown 
below. 


u 

cos  0 

sin#  sin  0 

sin#  cos  (j) 

V 

V 

= 

0 

cos  <j) 

-sixuj) 

w 

-sin# 

cos#  sin  0 

COS# COS  (j) 

_WC_ 

4.2  Preparing  Input  Data 

The  measured  airwake  data  file  which  is  used  as  input  to  AWSM  must  be  in  the  format 
defined  below: 


Column  1: 

Time 

Columns  2-4: 

um’vm’wm  f°r  the  top  anemometer  array 

Columns  5-7: 

um  >  vm  ’  wm  f°r  the  mid  anemometer  array 

Columns  8-10: 

um’vnvwm  f°r  the  low  anemometer  array 

Columns  11-12: 

V ship 5  ^  ship 

Columns  13-14 

Vref^ref 

Columns  15-16: 

Y stbd  > stbd 

Columns  17-18: 

V 

v  port’  1  port 

In  order  to  create  this  input  file,  a  number  of  utility  programs  (written  in  Fortran)  are 
available  for  use  on  Macintosh  computers. 

Program  MacTrans  is  first  used  to  obtain  data  in  imperial  units  from  the  data  file  generated 
using  program  MacShipRefine.  Appendix  C  contains  a  copy  of  file  AIRWAKE .  BLK,  which 
gives  the  channels  required  in  the  correct  order  for  use  with  utility  program  Convert  outlined 
below.  Before  running  MacTrans,  file  AIRWAKE .  BLK  should  be  renamed  to  TRANS  .  BLK. 
MacTrans  is  then  run  as  shown  below: 

[TRANS  version  date  12-FEB-91] 

GO  STRAIGHT  TO  EVENT  LOOP? 

I/P  FILENAME  =  20091158 

Wind  60  deg,  20  kn  -  Position  10 

File  20091158  (filtered) 

I/P  FILE  RECORDED  ON  17-FEB-93  AT  10:11:28 

INTEGN  INT  =  0.0000E+00;  RUN  CPU  TIME  =  3  MIN.  24.47  SEC. 

TIME  FROM  0 . 0000E+00  TO  8.9500E+01  IN  STEPS  OF  5.0000E-02 

*  PRC 

PRINTING  IN  COLUMNS  : 

BLKS 

L1,L2,L3 

<CR> 

IS  O/P  TO  TTY  REQRD  :  N 

[  PRC  Output,  for  this  run,  going  to  DSK:20091158 .COO  ] 

*  GOE 

**  RUNNING  ** 

*EXI 
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The  resulting  data  file,  20091158.  COO,1 2  is  then  used  as  input  to  program  AWSM. 

The  utility  program  Convert  may  be  used  to  change  the  MacTrans  column  format  into 
multiple  sequential  columns,  which  is  also  more  easily  read  in  by  other  programs.  Running 
program  Convert  results  in  file  2  0  0  9 1 1 5  8  .  COO  Ou  t .  The  operation  of  this  program  is 
similar  to  FFGConvert( Attitudes) 2  (see  Section  3.2). 

4.3  Running  AWSM 

When  program  Compat  creates  an  output  file  containing  calculated  motion  data,  the  file  is 
always  named  var .  dat.  When  multiple  runs  with  different  data  files  are  used,  it  is 
suggested  that  each  file  be  renamed  depending  on  the  run.3  For  the  example  used  throughout 
this  report,  the  var .  dat  file  is  renamed  2  0091158var .  dat. 

Program  AWSM  begins  by  flashing  an  FFG-7  picture  to  the  screen.  The  user  clicks  once  on 
this  picture  to  bring  a  window  to  the  screen  prompting  for  the  name  of  the  Compat  output  file 
(Fig.  12). 


Figure  12.  Window  Prompt  for  Ship  Motion  Data  File 


After  choosing  the  desired  motion  data  file,  the  user  is  then  prompted  for  the  measured 
anemometer  data  file  (Fig.  13). 


1  Program  MacTrans  has  an  automatic  file  naming  feature.  If  a  xxxxxxxx .  COL  file  already  exists  in  the 
same  folder,  the  new  file  will  be  named  xxxxxxxx .  COO  and  so  on  up  to  xxxxxxxx .  C9  9. 

2  Programs  Convert,  FFG  Convent  Attitudes),  and  FFGConvert(Heave)  are  basically  the  same  program  with 
minor  changes  included  in  data  manipulation. 

3  If  a  large  number  of  runs  are  to  be  made,  the  source  code  for  Compat,  namely  compat .  f ,  should  be  edited 
to  automatically  name  the  var .  dat  file. 
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Figure  13.  Window  Prompt  for  Anemometer  Data  File 


The  user  is  then  prompted  for  the  position  of  the  mobile  anemometer  file  for  that  particular 
run  (in  this  case  position  10),  as  shown  in  Fig.  14. 


Figure  14.  Window  Prompt  for  Anemometer  Position 


The  resulting  data  file,  in  this  case  2  0091158  Output,  contains  airwake  data  corrected  for 
motion  effects  and  may  be  compared  directly  with  the  uncorrected  data  in  file 
20091158  .COO  Out  (Section  4.2).  An  example  of  such  a  comparison  is  given  in  Section 
5.3. 

5.  RESULTS 

5.1  Motion  Analysis  Sensitivity  to  Data  Time  Length 

It  is  generally  accepted  that  analysis  of  ship  motion  requires  data  measured  over  a  long  time 
period.  During  the  airwake  trial,  it  was  impractical  to  record  data  at  each  of  the  anemometer 
positions  for  much  more  than  90  seconds.1  To  check  the  validity  of  the  assumptions  used  in 
the  motion  analysis  over  relatively  short  time  periods,  it  was  decided  to  examine  the  results 
using  the  same  data  file  for  varying  time  lengths.  The  file  20091347  was  chosen,  as  the  data 
were  recorded  for  118  s  (considerably  longer  than  the  usual  90  s)  and  the  magnitude  of  the 
ship  motion  for  this  case  was  relatively  large  (see  Appendix  D). 

Figures  15  and  16  show  comparisons  of  the  calculated  attitudes  and  velocities  respectively 
for  data  time  lengths  of  60,  75,  90, 105,  and  118  seconds. 


1  In  order  to  obtain  some  of  the  required  high  wind-over-deck  velocities  it  was  necessary  for  the  ship  to  steam 
at  full  speed.  Due  to  the  increase  in  fuel  consumption  at  this  speed,  the  data  were  recorded  at  each  of  the 
anemometer  positions  for  a  maximum  of  60  seconds  rather  than  the  usual  90  seconds. 
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As  can  be  seen  from  Fig.  15,  the  calculated  attitudes  give  similar  results  throughout  the  time 
periods  examined,  indicating  that  the  method  for  determining  the  attitudes  is  relatively 
insensitive  for  data  recorded  for  60  s  or  more.  Examination  of  Fig.  16  shows  that  the  vertical 
velocity  component  gives  similar  results  for  the  different  time  periods,  but  the  longitudinal 
component  gives  slight  differences  and  the  lateral  component  varies  for  each  data  length. 
These  differences  suggest  that  some  of  the  assumptions  used  for  determining  the  velocities 
may  not  be  good.  Some  indication  of  the  confidence  level  of  the  estimates  is  given  by  the 
Cramer- Rao  bounds  for  each  of  the  parameters.  Table  6  gives  the  Cramer-Rao  bounds  for 
the  parameters  estimated  for  each  of  the  different  data  time  periods  and  Fig.  17  shows  a 
comparison  of  the  Cramer-Rao  bounds  for  the  accelerometer  offsets. 


Table  6.  Cramer-Rao  Bounds  for  Various  Data  Time  Periods 


Parameter 

Cramer-Rao  Bounds 

Name 

No. 

60s 

75  s 

90s 

105  s 

118s 

Pb 

4 

0.000031 

0 . 000034 

0.000025 

0.000020 

0.000016 

% 

5 

0.000010 

0 . 000007 

0 . 000005 

0 . 000004 

0 . 000003 

rb 

6 

0.000001 

0.000001 

0.000001 

0.000000 

0 . 000000 

K 

28 

0.001069 

0.001467 

0.001302 

0.001225 

0.001092 

Vo 

29 

0.000349 

0.000311 

0.000268 

0.000241 

0.000223 

axb 

1 

0.000638 

0.000391 

0.000364 

0.000225 

0.000244 

ayb 

2 

0.001325 

0.001244 

0.001361 

0.001424 

0.000863 

azb 

3 

0.000340 

0.000218 

0.000133 

0 . 000086 

0 . 000063 

Uo 

25 

0.021544 

0.018110 

0.016181 

0.013652 

0.013135 

Vo 

26 

0.045881 

0.053726 

0.070686 

0.086244 

0.058998 

27 

0.010628 

0.008361 

0.006167 

0.004743 

0.003850 

38 

0.139626 

0.136450 

0.120662 

0.106561 

0.098713 
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Vertical  Velocity  (ft/s)  Lateral  Velocity  (ft/s)  Longitudinal  Velocity  (ft/s) 


Figure  16.  Comparison  of  Calculated  Ship  Velocities  for  Data  of  Various  Time  Periods 
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Figure  17.  Comparison  of  Cramer-Rao  Bounds  for  Accelerometer  Offsets 

It  can  be  seen  from  Fig.  17  that  the  C-R  bounds  for  the  estimates  of  the  longitudinal  and 
vertical  accelerometer  offsets,  aXb  and  aZb  (which,  in  effect,  determine  the  calculated 

velocities),  steadily  decrease  with  an  increase  in  the  length  of  recorded  data.  This  suggests 
that  the  basic  assumptions  are  correct  and  the  longer  the  length  of  data  recorded,  the  higher 
the  confidence  level  for  the  estimated  offsets.  However,  the  C-R  bound  for  the  lateral 
accelerometer  offset,  av  ,  is  much  larger  than  that  for  the  other  two  offsets  and  does  not 

show  a  consistently  decreasing  trend  with  time.  This  suggests  that  the  assumption  that  the 
average  lateral  velocity  is  zero  may  be  in  error,  perhaps  due  to  ocean  currents.  To  improve 
the  accuracy  of  this  technique  in  determining  the  lateral  velocity  component,  the  actual  lateral 
velocity  or  displacement  could  be  measured,  e.g.  by  using  a  satellite  navigation  system. 

5.2  Determining  Worst  Case  Ship  Motion 

Before  proceeding  with  determining  ship  motion  and  removing  these  effects  from  all  the  data 
files  (a  total  of  146  data  runs),  it  was  decided  to  use  Compat  to  find  the  approximate 
magnitude  of  the  motion  for  each  of  the  wind-over-deck  conditions.  Once  the  worst  case  (i.e. 
largest  linear  velocities  and  angular  rates)  was  found,  the  airwake  without  ship  motion  would 
be  determined  for  that  complete  case.  Table  7  gives  a  list  of  the  files  which  were  initially 
processed. 

The  order  in  which  the  files  in  Table  7  were  processed  was  largely  historical,  as  discussed 
next.  By  examining  the  raw  data,  a  number  of  possibilities  presented  themselves  as  the  worst 
case.  The  initial  examinations  were  done  for  the  mobile  mast  positioned  over  the  bullseye 
(position  8).  At  this  stage,  the  180°  10  kn  wind-over-deck  cases  appeared  to  be  giving  the 
greatest  magnitudes  of  ship  motion,  so  these  files  were  concentrated  upon  to  finalise  the 
technique  used  to  determine  ship  motion.  Once  the  technique  was  finalised,  a  representative 
file  from  each  of  the  wind-over-deck  conditions  which  had  not  yet  been  examined  was 
processed  for  a  deck  position  resulting  in  the  greatest  effect  on  the  anemometers  (position  10 
or,  if  that  was  not  available,  position  5).  It  should  be  noted  that  for  the  90°  20  kn  case  (where 
HF  radio  corrupted  many  of  the  digital  channels),  it  was  found  that  the  files  for  positions  5  1 


1  Position  10  was  not  measured. 
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and  1  (21091403  and  21091409  respectively)  also  had  errors  in  the  yaw  rate  measurements,1 
which  resulted  in  incorrect  ship  motion  results. 


Table  7.  Initial  Selected  Wind-Over-Deck  Conditions 


Wind-Over-Deck 

Condition 

File  Name 

Mobile  Mast 
Position 

Ship  Stationary 

16091506 

8 

00°  1 0  kn 

19091053 

10 

00°  20  kn 

19091624 

8 

00°  35  kn 

19091714 

5 

30°  1 0  kn 

20091424 

8 

30°  20  kn 

18091428 

10 

30°  35  kn 

21091243 

8 

60°  10  kn 

18091529 

10 

60°  20  kn 

20091158 

10 

90°  10  kn 

18091626 

8 

90°  20  kn 

21091424 

11 

135°  10  kn 

20091350 

5 

135°  10  kn 

20091315 

10 

180°  10  kn 

20090930 

6 

180°  10 kn 

20090956 

8 

180°  10 kn 

20091000 

9 

180°  10 kn 

20090920 

13 

An  Excel  (©Microsoft)  document  (Stats  Template)  was  set  up  to  find  the  maximum, 
minimum,  mean,  and  standard  deviation  of  the  ship  motion  variables.  The  output  files  from 
Compat  are  pasted  directly  into  the  first  14  columns  of  the  Excel  document  with  the  statistics 
being  given  in  the  next  14  columns.  After  examination  of  all  the  files,  it  was  determined  that 
the  worst  ship  motion  occurred  for  the  180°  10  kn  and  the  135°  10  kn  cases.  At  this  stage  it 
was  decided  to  concentrate  on  the  135°  10  kn  case,  as  this  seemed  to  have  more  reasonable 
motion  data  than  the  180°  10  kn  case.2  The  files  for  the  135°  10  kn  case  are  shown  in 
Table  8. 

Statistics  for  the  135  deg  10  kn  case  are  presented  in  Appendix  D  with  statistics  for  the 
remaining  files  in  Table  7  given  in  Appendix  E. 

5.3  Comparing  Measured  with  Corrected  Airwake  Data 

Having  determined  the  ship  motion  for  the  entire  135°  10  kn  case,  program  AWSM  was  used 
to  determine  the  airwake.  Table  9  gives  the  statistics  of  the  top  anemometer  array  for 
position  5,  comparing  measured  data  with  data  corrected  for  motion  effects,  while  Figure  18 
shows  the  time  histories  of  the  lateral  and  vertical  components.  The  corrections  for  the 
longitudinal  velocity  component  for  this  case  were  so  small  as  to  be  negligible. 


1  The  errors  in  the  yaw  rate  measurements  consisted  of  a  large  step  part  way  through  the  time  history,  which 
could  possibly  be  corrected  after  further  analysis. 

2  The  vertical  displacements,  Z,  for  the  180°  10  kn  case  appear  to  be  abnormally  large  when  compared  to  all 
the  other  cases.  This  may  however  be  true  because  the  ship  was  almost  stationary  to  achieve  10  kn  wind- 
over-deck  from  behind,  and  thus  the  ship  was  not  cutting  through  the  waves  but  rather  wallowing  on  them. 
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Examination  of  Fig.  18  shows  that,  at  a  given  point  in  the  time  history,  the  magnitude  of  the 
airwake  velocity  when  corrected  for  ship  motion  differs  markedly  from  that  when 
uncorrected,  although  the  mean  velocity  is  almost  identical. 

Table  8.  Files  for  135°  10  kn  Case 


Filename 

Mast 

Position 

20091337 

1 

20091340 

2 

20091343 

3 

20091347 

4 

20091350 

5 

20091330 

6 

20091334 

7 

20091307 

8 

20091311 

9 

20091315 

10 

20091326 

11 

20091323 

12 

20091319 

13 

Table  9.  Statistics  for  Top  Anemometer  Array  at  Position  5  for  135°  10  kn  Case 


Statistics 

um  (ft/s) 

u  (ft/s) 

vm  (ft/s) 

v  (ft/s) 

wm  (ft/s) 

w  (ft/s) 

Minimum 

5.37 

5.77 

-18.62 

-21.10 

-6.78 

-8.06 

Maximum 

12.91 

13.87 

-7.91 

-6.97 

0.00 

1.51 

Mean 

10.65 

10.67 

-11.56 

-11.52 

-3.75 

-3.85 

Std  Deviation 

1.28 

1.54 

2.07 

2.52 

0.97 

1.59 

The  maximum  corrections  to  the  lateral  velocity  component  are  approximately  3.3  ft/s  for 
ship  linear  velocity  ( AVship )  and  2.0  ft/s  ( pAziat )  due  to  rolling.  The  time  at  which  the 

correction  for  ship  linear  velocity  is  maximum  is  generally  when  the  correction  due  to  rolling 
is  minimum,  and  vice  versa.  Even  though  these  corrections  have  been  applied,  Fig.  18 
indicates  a  signal  having  approximately  the  same  period  as  the  ship  motion.  Preliminary 
study  of  a  number  of  cases  suggests  that  this  is  due  to  large  eddy  turbulence  in  the  earth’s 
boundary  layer,  rather  than  vortices  being  shed  from  the  sharp  edges  of  the  ship. 
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Figure  18.  Comparison  of  Measured  Data  with  Data  Corrected  for  Motion  Effects 
6.  CONCLUDING  REMARKS 

A  procedure  has  been  presented  here  for  determining  ship  motion  from  limited  measured  data 
and  removing  this  effect  from  measured  airwake  data.  The  same  technique  for  determining 
ship  motion  may  be  used  in  future  trials,  with  a  notable  difference  being  that  if  attitude  is 
measured  directly,  then  it  is  not  necessary  to  use  the  mean  values  derived  from 
accelerometers.  The  accuracy  of  the  calculated  ship  velocity  may  be  improved  by 
lengthening  the  data  recording  time  and  by  including  extra  velocity  and  displacement 
information  from  a  satellite  navigation  system. 

It  has  been  shown  that  the  mean  of  the  measured  airwake  is  approximately  the  same  as  the 
mean  of  the  airwake  corrected  for  motion  effects.  With  hindsight,  this  is  an  inevitable  result, 
as  the  mean  angular  rates  must  approach  zero  and.  the  ship  is  assumed  to  have  mean  linear 
velocities  approaching  zero  in  the  motion  analysis.  For  the  applications  of  using  the 
measured  data  as  a  mean  flow  data  base  for  the  helicopter/ship  model,  or  for  comparison  with 
wind  tunnel  results,  the  procedure  for  removing  motion  effects  from  the  airwake  is  therefore 
not  warranted.  However,  it  should  be  noted  that  the  technique  for  determining  ship  motion 
developed  here  has  also  been  applied  to  determining  the  motion  of  a  ‘Perth’  class  DDG 
destroyer  and  for  calculating  the  position  of  the  flight  deck  of  HMAS  Jervis  Bay  during  a 
Black  Hawk  First  of  Class  Flight  Trial. 
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APPENDIX  A  -  INPUT  FILES  COMDAT.  002  AND  COMDAT.  004 


Parameter 

Definition 

Parameter 

Definition 

NPM 

No.  of  parameters 

ITR 

No.  of  iterations 

MU 

No.  of  control  inputs 
(=9;  ax,ay,az,p,q,r,p,q,r) 

ITO 

No.  of  iterations  with  a  priori 
values  applied 

NS 

No.  of  states 

(=9;  U,V,W,<j),e,V,X,Y,Z ) 

NAVG 

No.  of  points  used  for  least 
squares  fit  of  initial  data  point 

MZ 

No.  of  measurements  (=9; 

^ m ’ »  W/n >  fim >  ■> 

DTT 

Data  time  interval 

ISK 

NN 

No.  of  lines  skipped  in  data  file 

No.  of  lines  in  data  file 

FINK 

Program  variable 

COMDAT.002 


NEM 

MJ 

IB 

IE  I3C 

NSJ 

44 

9 

9 

9  2 

2500 

DIT 

EUK 

0.05 

0.001 

SIKIE3  RETS  SCflUS 

U 

1 

1 

v 

2 

1 

w 

3 

1 

ehi 

4 

1 

HEIR 

5 

1 

ESI 

6 

1 

X 

7 

1 

Y 

8 

1 

Z 

9 

1 

MESS 

KEEM  SEGUG 

I.C.  D1 

U 

1 

1 

0  0.000001 

v 

2 

1 

0  0.000001 

w 

3 

1 

0  0.000001 

hh 

4 

1 

0  100000.0 

iHEim 

5 

1 

0  100000.0 

psr 

6 

1 

0  100000.0 

X 

7 

1 

0  0.000001 

Y 

8 

1 

0  0.000001 

Z 

9 

1 

0  0.000001 

EHR2TE 

KEEP  SIMUS 

A  HUCRI 

VEiair 

BAX 

1 

0 

0.00 

99999.0 

BAY 

2 

0 

0.00 

99999.0 

BAZ 

3 

0 

0.00 

99999.0 

BP 

4 

1  -0.000000 

99999.0 

B2 

5 

1  -0.000000 

99999.0 

m 

6 

1 

0.000000 

99999.0 

IAX 

7 

0 

0.00 

99999.0 

IAY 

8 

0 

0.00 

99999.0 

IAZ 

9 

0 

0.00 

99999.0 

IP 

10 

0 

0.00 

99999.0 

IQ 

11 

0 

0.00 

99999.0 

IR 

32 

0 

0.00 

99999.0 

BU 

13 

0 

0.00 

99999.0 

W 

14 

0 

0.00 

99999.0 

EW 

15 

0 

0.00 

99999.0 

bhu 

16 

0 

0.00 

99999.0 

EfflEEA 

17 

0 

0.00 

99999.0 

BPSI 

18 

0 

0.00 

99999.0 

m 

19 

0 

0.00 

99999.0 

IV 

20 

0 

0.00 

99999.0 

IW 

21 

0 

0.00 

99999.0 

lhh 

22 

0 

0.00 

99999.0 

TflWRm 

23 

0 

0.00 

99999.0 

IPSI 

24 

0 

0.000000 

99999.0 

rro  m 

2  10 


Al 


uo 

25 

0 

0.000000 

99999.0 

VO 

26 

0 

0.00 

99999.0 

wo 

27 

0 

0.00 

99999.0 

HUO 

28 

1  -0.000000 

99999.0 

wo 

29 

1  -0.000000 

99999.0 

PStO 

30 

0 

0.00 

99999.0 

BX33 

31 

0 

-87.7 

99999.0 

EftCG 

32 

0 

4.33 

99999.0 

BZCG 

33 

0 

-14.9 

99999.0 

HUHT 

34 

0 

0.00 

99999.0 

H3HT 

35 

0 

0.00 

99999.0 

XD 

36 

0 

0.00 

99999.0 

YO 

37 

0 

0.00 

99999.0 

ZO 

38 

0 

0.00 

99999.0 

BX 

39 

0 

0.00 

99999.0 

BY 

40 

0 

0.00 

99999.0 

EZ 

41 

0 

0.00 

99999.0 

m 

42 

0 

0.00 

99999.0 

w 

43 

0 

0.00 

99999.0 

w ; 

44 

0 

0.00 

99999.0 

COMDAT.004 

NHM 

MJ 

N3 

mz  iac  m 

44 

9 

9 

9 

2  2500 

DIT 

EMC 

0.05 

0.001 

SIKHS  REFS  SIKTOB 

U 

1 

1 

V 

2 

1 

W 

3 

1 

mi 

4 

1 

THEIA 

5 

1 

PSI 

6 

1 

X 

7 

1 

Y 

8 

1 

Z 

9 

1 

Mas 

KEFM  SIKEUS 

I.C. 

D1 

u 

1 

1 

0  100000.0 

V 

2 

1 

0  100000.0 

w 

3 

1 

0  100000.0 

mr 

4 

1 

0  0.000001 

HJEIK 

5 

1 

0  0.000001 

PSI 

6 

1 

0  0.000001 

X 

7 

1 

0  100000.0 

Y 

8 

1 

0  100000.0 

Z 

9 

1 

0  100000.0 

PARSM3  KEEP  SIKIU3 

A  PRICRI 

WEIGHT 

BSC 

1 

1 

0.109218 

99999.0 

BHY 

2 

1 

1.247285 

99999.0 

mz 

3 

1  -0.120199 

99999.0 

BP 

4 

0  -0.020238 

99999.0 

B2 

5 

0 

0.000979 

99999.0 

HI 

6 

0 

0.000044 

99999.0 

IAX 

7 

0 

0.00 

99999.0 

iay 

8 

0 

0.00 

99999.0 

IAZ 

9 

0 

0.00 

99999.0 

BP 

10 

0 

0.00 

99999.0 

IQ 

11 

0 

0.00 

99999.0 

LR 

12 

0 

0.000000 

99999.0 

HJ 

13 

0 

0.00 

99999.0 

BJ 

14 

0 

0.00 

99999.0 

EW 

15 

0 

0.00 

99999.0 

BEHI 

16 

0 

0.00 

99999.0 

BIHE3K 

17 

0 

0.00 

99999.0 

BPSI 

18 

0 

0.00 

99999.0 

HJ 

19 

0 

0.00 

99999.0 

IV 

20 

0 

0.00 

99999.0 

IW 

21 

0 

0.00 

99999.0 

BEHI 

22 

0 

0.00 

99999.0 

HtHEIA 

23 

0 

0.00 

99999.0 

IRSI 

24 

0 

0.00 

99999.0 

UO 

25 

1 

7.8292 

99999.0 

VO 

26 

1 

0.00 

99999.0 

WO 

27 

1 

0.00 

99999.0 

MIO 

28 

0  -0.006754 

99999.0 

N®G 

10 


A2 


TOO 

29 

0  -0.008442 

99999.0 

RSIO 

30 

0 

0.00 

99999.0 

BXG3 

31 

0 

-87.7 

99999.0 

BH33 

32 

0 

4.33 

99999.0 

ECG 

33 

0 

-14.9 

99999.0 

Bomr 

34 

0 

0.00 

99999.0 

H3HI 

35 

0 

0.00 

99999.0 

XD 

36 

0 

0.00 

99999.0 

YO 

37 

0 

0.00 

99999.0 

ZO 

38 

1 

-2.317998 

99999.0 

EK 

39 

0 

0.00 

99999.0 

W 

40 

0 

0.00 

99999.0 

EZ 

41 

0 

0.00 

99999.0 

m 

42 

0 

1.00 

99999.0 

w 

43 

0 

1.00 

99999.0 

w. 

44 

0 

1.00 

99999.0 

APPENDIX  B  -  PROBLEM  SPECIFIC  SUBROUTINES  AND  PROGRAM  LINK  FILE 


Link 

xlf  -g  carpat.f  cxmsub.f  deriv.f  init.f  re^n.f  rk4new.f  -o  CCrrpatH 


Deriv.f 

SHEUIEKE  UREVS(  F,  NK,  PARAM,  UJ,  iffiR,  N3TS  ) 

C 

C  SUHCUHNE  KR  INIHRAffENS  AIRCRAFT  KINEMACCrC  FQFKnCFB 
C  KR  6DCF  (XMPAT1BILJIY  CHECKIIC 
C 

REAL  AX,  AY,  AZ,  CUAR4,  CUAR5,  F(  *  ),  P,  EARAM(  *  ) 

REAL  Q,  R,  SUAR4,  S®R5,  UJ(  *  ) ,  VAR(  FEES,  *  ) 

C 

cAMA  (RAV=9. 80665  !  m/s/s 

EARAM3IER{  (RAV=32.17405  )  !  ft/s/s 

c 

cAMA  Accelerations  and  rates  are  read  in  as  cmtrol  inputs  fran  a  (feta 
cAMA  file  created  hy  'FPNrP.FEG'  virLch  nanipulates  feta  fran  a  'TRANS' 
c®R  colum  print 
c 

cAMA  Included  angular  accelerations  derived  fran  filtered  rates 
C 

P  =  (  1.0  +  EARAM(  10  )  )*UU(  4  )  +  EARAM(  4  ) 

Q  =  (  1.0  +  PARAM(  11  )  )*UJ(  5  )  +  EKRAM(  5  ) 

R  =  (  1.0  +  EARAM(  12  )  )*UCJ(  6  )  +  PARAM(  6  ) 

pdot  =  uu(7) 
qlot  =  uu(8) 
rfet  =  uu(9) 
c 
C 

SUAR4  =  SIN(  VAR(  4,  INK  )  ) 

CVAR4  =  CDS(  VAR(  4,  NK  )  ) 

SUAR5  =  SIN(  IffiR(  5,  IK  )  ) 

CVAR5  =  CC6(  VAR(  5,  KK  )  ) 
svar6  =  sin(  var(  6,  rik  )  ) 
cvar6  =  cos(  var(  6,  nk  )  ) 
c 

Ax  =  UU(  1  )  +  (  Q*Q  +  R*R  )  *PARAM(  31  )  -  (P*Q-rfet)  *PARAY[(  32  ) 

$  -  (P*RKj±>t)*EARAM(  33  )  -  GRAV*S7AR5 

Ay  =  UU(  2  )  -  (P*Qfrdot)  *PARAM(  31  )  +  (  P*P  +  R*R  )  *EARAM(  32  ) 

$  -  (Q*R-pcbt)  *PARAM(  33  )  +  (RAV*OffiR5*S\ffiR4 

Az  =  U(J(  3  )  -  (P*R-feot)  *PARAM(  31  )  -  (Q*Rnpdot)  *PARAM (  32  ) 

$  +  (  P*P  +  Q*Q  )*PARAM(  33  )  +  GRA5/*ENffiR5*CVAR4 

C 

AX  =  (  1.0  +  EARAM(  7  )  )*Ax  +  ERRAM(  1  ) 

AY  =  (  1.0  +  PARAM(  8  )  )*Ay  +  EHRAM(  2  ) 

AZ  =  (  1.0  +  PARAM(  9  )  )*Az  +  PARAM(  3  ) 

C 

C  CHLOULAIIE  EEREVAEIVES 
C 

cAMA  Eerivative  agnations  (Ahich  are  integrated  to  give  cteervaticns 
cAKA  Elgjaticns  are  lEot,VDot,VCDt,EhiDDt,EhetaD3t,IfeinDt,XDDt,YDDt,ZDDt 
C 

c  vel  =  sqrt(var(l,rik)**2  +  var(2,nk)**2 
c  $  +  var(3,nk)**2) 

F(  1  )  =  AX  -  Q*\ffiR(  3,  NK  )  +  R*UAR(  2,  NK  )  !  U  DOT 

F(  2  )  =  AY  -  R*VAR(  1,  KK  )  +  P*\ffiR(  3,  NK  )  !  V  EOT 

F(  3  )  =  AZ  -  P*VAR(  2,  NK  )  +  Q*VAR(  1,  NK  )  !  W  DOT 

F(  6  )  =  (  Q*SVAR4  +  R*CVAR4  )/OZAR5  !  PSI  EOT 

F(  4  )  =  P  +  F(  6  )*SVAR5  !  HU  EOT 

F(  5  )  =  Q*CWAR4  -  R*S\fflR4  !  EHEEA  DOT 

cAMA  f  (  7  )  =  vel*cvar5*cvar6*param(42)  !  X  Dot 

cAMA  f  (  8  )  =  vel*cvar5*svar6*param(43)  !  Y  Dot 

cAMA  f  (  9  )  =  -vel*svar5*paran(44)  !  Z  Eot 

f(  7  )  =  (var(l,rik)  *cvar5*cvar6 
$  +  var(2,nk)*(svar4*svar5*cvar6-cvar4*svar6) 

$  +  var(3,nk)  *  (cvar4*svar5*cvar6+svar4*svar6) )  *paran(42) 

f(  8  )  =  (var(l,nk)*cvar5*svar6 
$  +  var(2,nk)*(svar4*svar5*svar6+cvar4*cvar6) 

$  +  var  (3,nk)  *  (cvar4*svar5*svar6-svar4*cvar6) )  *param(43) 


B1 


f(  9  )  =  (-var(l,nk)*svar5  +  var(2,nk)*cvar5*svar4 
$  +  %^(3,nk)*cvar5*cvar4)*param(44)  !  Z  Dot 

RETCFN 

END 


Init.f 


SUEKUnNE  EARSMJEHTf  EARAM,  CUT,  IC_M,  MG,  Z,  ZO,  MZ,  MZTS  ) 

* 

*  Subroutine  to  set  initial  axditicn  parameters . 

* 

*  Modified  h/:  M.L.  ’Dimer,  SDG,  ARL 

*  Date:  4/6/ '91 

•k 

k _ 

c 

C  SUERCUITEE  ID  ST  INTITAL  (XEDIITCES  ECR  PARSMEIEEiS 
C 

C  FCEICWTEC  ALICEE  ECR  EHHER  HENHElCmCN  CF  MITTAL 
C  CCEDIITCES,  CR  ECR  UEEKITN3  BASED  CN  HENITFTED  BIAS 
C  AED  S3LE  F7CICR  PARAMETERS  -  ASSLMEN3  SIKHS  ARE 
CU(0),  V(0),  W(0),  EHI(O) ,  THEIK(O),  PSI(O) 

C 

INIHHR  1C_M(  *  ) ,  I®M3,  EG,  EEIS 
C 

REAL  AV3,  UN,  DIT,  NUM,  PARAM(  *  ) ,  SU4TTEE,  SJ4_TIME_® 

REAL  THE,  TTEE_Z0(  20  ),  Z(  EEIS,  *  ),  Z0(  *  ) 

C 

AU3  =  Fimr(  EKV3  ) 

DO  I  =  1,  MZ 
ZQ(  I  )  =  0.0 

TIME_Z0(  I  )  =  0.0 
ENID 

StMJITEE  =  0.0 
Sfe_TTME_S2  =0.0 
C 

DO  J  =  1,  NAV3 

TIME  =  DIT*ElDKr(  J  -  1  ) 

StMTTME  =  StMTTME  +  TIME 
SLMJTME-92  =  SCM_TTME_S2  +  TIME*TIME 
DO  I  =  1,  IE 

3F(  IC_M(  I  )  -NE.0  )  THEN 
Z0(  I  )  =  Z0(  I  )  +  Z(  I,  J  ) 

TIME_Z0(  I  )  =  TTME_Z0(  I  )  +  TIME*Z(  I,  J  ) 

ENDIF 

ENCDO 

EEDDO 

UN  =  A\G*  StM_TTME__S2  -  SCMJITME*SLM_TTME 
C 

DO  I  =  1,  EE 

IF(  IC_M(  I  ).NE.O  )  THEN 
NUM  =  A\G*TTME_Z0(  I  )  -  Z0(  I  )*SM_TME 
Z0(  I  )  =  (  Z0(  I  )  -  StMJTEETDM/EEN  )/AV3 
ENDIF 
ENID 
C 

C  Initial  oxditicns  far  u, phi, theta, psi  rpfeted  since  (feta  measure!, 

C  tut  IQs  far  v,w  estimated  since  the/  are  unkrxan. 

C 

IF(  ICM(  1  ).EE.O  )  THEN 

PARSM(  25  )  =  (  Z0(  1  )  -  EARAM(  13  )  )/(  1.0  +  PAR?M(19  )  )  !  U0 
EEDIF 

IF(  IC_M(  2  ).NE.O  )  THEN 

EARfiM(  26  )  =  (  Z0(  2  )  -  PARAM(  14  )  )/(  1.0  +  PARAM(  20  )  )  !  VO 
EEDIF 

IF(  IC_M(  3  ).NE.O  )  THEN 

PARAMf  27  )  =  (  Z0(  3  )  -  PAR?M(  15  )  )/(  1.0  +  EARSM(  21  )  )  i  W0 
EEDIF 

IF(  IC_M(  4  )  .EE.O  )  THEN 

EARAM(  28  )  =  (  Z0(  4  )  -  PARAM(  16  )  )/{  1.0  +  PARAMf  22  )  )  !  TtetaO 
EEDIF 

IF(  IC_M(  5  )  .EE.O  )  THEN 


B2 


BKRPM(  29  )  =  (  Z0(  5  )  -  EARAM(  17  )  )/(  1.0  +  EARAM(  23  )  )  !  flliO 
EN3EF 

IF(  ICM(  6  ).NE.O  )  THEN 

ERF@M(  30  )  =  (  Z0(  6  )  -  EARAM(  18  )  )/(  1.0  +  EARAM(  24  )  )  !  EsiO 
BSDIF 
C 

if (  ic_m(  7  ).ne.O  )  then 
parsm(  36  )  =  (  z0(  7  )  -  param(  39  )  )  !  XD 

axlif 

if(  ic_m(  8  )  .ne.O  )  thai 
paran(  37  )  =  (  z0(  8  )  -  param(  40  )  )  !  YO 

axlif 

if  (  icm(  9  )  .ne.O  )  thai 
param(  38  )  =  (  zO(  9  )  -  paran(  41  )  )  !  ZO 

svdif 
C 

RE3UFN 

END 

C 


SCMBOITNE  SMTEjmrrt  VAR,  IK,  EARAM,  IBIS  ) 

* 

*  Subroutine  to  set  initial  cxnditicns  far  state  equations. 

★ 

*  Modified  by:  M.L.  Turner,  STS,  ARL 

*  Date:  16/5/ '91 

* 

* _ 

C 

C  The  following  allows  far  either  idmtificaticn  of  initial  ccnditicns 

C  or  far  rpcbting  based  cn  idsitified  bias  and  scale  factor  parameters 

C 

REAL  EAR*[(  *  ),  VAR(  IBIS,  *  ) 

C 

C  PARAM(  25  )  =  (  Z0(  1  )  -  EARAM(  13  )  )/(  1.0  +  PARAM(  19  )  )  !  IX) 

C  PARAM(  26  )  =  (  Z0(  2  )  -  PAR£M(  14  )  )/(  1.0  +  PARAM(  20  )  )  !  V0 

C  PARAM(  27  )  =  (  Z0(  3  )  -  EARAM(  15  )  )/(  1.0  +  PARAM{  21  )  )  !  WO 

C  PARAM(  28  )  =  (  Z0(  4  )  -  EARAM(  16  )  )/(  1.0  +  PARAM(  22  )  )  !  ThetaO 

C  PARAM(  29  )  =  (  Z0(  5  )  -  EffiRAM(  17  )  )/(  1.0  +  PARAM(  23  )  )  !  RiiO 

C  PARAM(  30  )  =  (  Z0(  6  )  -  EARAM(  18  )  )/(  1.0  +  PARAM(  24  )  )  !  E£i0 

C 

VAR(  1,  IK  )  =  PARAM(  25  )  !  Estimates  uO 

VAR(  2,  IK  )  =  EARAM(  26  )  !  Estimates  vO 

VAR(  3,  IK  )  =  PARAM(  27  )  !  Estimates  wO 

VAR(  4,  NEC  )  =  PARAM  (  28  )  !  Estimates  phiO 

VAR(  5,  IK  )  =  EARAM(  29  )  !  Estimates  thetaO 

C  VAR(  6,  IK  )  =  PARAM  (  30  )  !  Estimates  psiO 

C 

AffiR(  6,  IK  )  =  0.0  !  PSIO  =  zero  by  definiticn 

c 

var(  7,  nk  )  =  0.0  !  xO  =  zero  by  definiticn 

var(  8,  nk  )  =  0.0  !  yO  =  zero  by  definiticn 

var(  9,  nk  )  =  param(  38  )  !  Estimates  zO 

C 

KE3URT 

END 

Resp.f 

SUHBOITNE  RESP(  AXRAW,  AYRAW,  P,  Q,  R,  Y,  IK,  PARAM,  UJ, 

$  VAR,  IBIS  ) 

C 

cAMA  Cbservaticns  vhich  are  catparedwith  Z  (maasunsrsits)  in  'HPl.Plot1 
cAMA 

cAMA  Ebdified  Ifcv  90  to  estimate  oartre  of  moticn  of  ship  and  allow 

cAMA  far  paxblum  attitude  instnrtHits . AE® 

C 

REAL  AXRAW,  AYRAW,  EelEtli,  DalTht,  PARAM(  *  ) ,  P,  Q,  R 
REAL  TXJ(  *  ),  VAR(  IBIS,  *  ) ,  Y(  *  ) 

C 

QOAGRAV  /  32.17405 
C 

O - 

C  UST  Y(l)  ,Y(2)  . .  .AS  EUCITCtB  CF  VAR(l) ,  VAR(2) 

C  EARAM(l) ,  EARAM (2)  ...  AND  UU(1)  ,UQ(2) . . . 


B3 


c - 

c 

c  VEL  =  S3*T(  VAR(  1,  NK  )**2  +  VSR(  2,  MC  )**2 

Cl  +  VAR(  3,  NK  )**2  ) 

C 

am  Correct  far  aoceleroreter  offset  -  PAR£M(31,32,33)  are  x,y,z  offsets 
C 

P  =  (  1.0  +  PARfflK  10  )  )*UJ(  4  )  +  4  ) 

Q  =  (  1.0  +  PAR£M(  11  )  )*TU(  5  )  +  EAR£M(  5  ) 

r  =  (  l.o  +  par?m(  12  )  )*uu(  6  )  +  mm(  6  ) 

AxRa w  =  uu(  1  )  +  (  Q*Q  +  R*R  )  *EftRAM(  31  ) 

$  -  (  P*Q  )  *P£F5W(  32  )  -  (  P*R  )  *PARfM(  33  ) 

$  -  OtAV*SIN(  VfR(  5,  UK  )  ) 

A>«aw  =  uu(  2  )  -  (  P*Q  )*PARSM(  31  ) 

$  +  (  P*P  +  R*R  )*mm(  32  )  -  (  Q*R  )*EJ*SM(  33  ) 

$  +  GtAV*CDS(  VSR(  5,  MC  )  )*SIN{  \ffiR(  4,  MC  )  ) 

c 

c  fimyi-  fca:  padfLun  attitude  instruimts  (pitch  and  roll) 
c 

EelTht  =  EHR?M(  34  )*AxRaw 

Ifellhi  =  PAR?M(  35  )*Ai*aw 

C 

Y(  1  )  =  (  1.0  +  PARfiMt  19  )  )*VSR(  1,  MC  )  +  PRR£M(  13  )  !  u  (ship  ^eed) 

Y(  2  )  =  (  1.0  +  PARSM(  20  )  )*VRR(  2,  MC  )  +  EftRSM(  14  )  !  v 

Y(  3  )  =  (  1.0  +  PAE@M(  21  )  )*VAR(  3,  MC  )  +  EM®M(  15  )  !  w 

am  Y(  4  )  =  (  1.0  +  EARAM(  22  )  )*VAR(  4,  MC  )  +  PAR£M(  16  )  +  tellhi  !  phi 

cSMA  Y(  5  )  =  (  1.0  +  PARSM(  23  )  )*VAR(  5,  MC  )  +  EAE®M(  17  )  +  DeMht  !  theta 

y(  4  )  =  (  1.0  +  param(  22  )  )*var(  4,  nk  )  +  paran(  16  )  !  phi 

y(  5  )  =  (  1.0  +  param(  23  )  )*var(  5,  nk  )  +  param(  17  )  !  theta 

Y(  6  )  =  (  1.0  +  PARSM(  24  )  )*VAR(  6,  MC  )  +  PARSM(  18  )  !  psi 

am  y(  7  )  =  var(7,rik)  +  paran(39)  !  X 

am  y(  8  )  =  var(8,rik)  +  paran(40)  !  Y 

y(  7  )  =  o.o 
y(  8  )  =  0.0 

y(  9  )  =  var(9,nk)  +  paran(41)  !  Z 
C 

REIUHSI 

BSD 


B4 


APPENDIX  C  -  MACTRANS  ‘.BLK’  FILES 
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23465  51  52  00000000000000000 

0  0  0  0  0  0 
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L2 

47  53  54  31  32  16  17  18  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

0  0  0  0  0  0 

L3 

19  00000000000000000000000 

0  0  0  0  0  0 
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APPENDIX  D  -  STATISTICS  FOR  135°  10  KN  RELATIVE  WIND 
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APPENDIX  E  -  SELECTED  SHIP  MOTION  STATISTICS 
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